BCS-BEC Crossover in Atomic Fermi Gases with a Narrow Resonance 
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We determine the effects on the BCS-BEC crossover of the energy dependence of the effective 
two-body interaction, which at low energies is determined by the effective range. To describe 
interactions with an effective range of either sign, we consider a single-channel model with a two- 
body interaction having an attractive square well and a repulsive square barrier. We investigate the 
two-body scattering properties of the model, and then solve the Eagles-Leggett equations for the 
zero temperature crossover, determining the momentum dependent gap and the chemical potential 
self-consistently. From this we investigate the dependence of the crossover on the effective range of 
the interaction. 
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I. INTRODUCTION 

Recently there has been remarkable progress in realiz- 
ing, in ultracold atomic Fermi gases, the crossover from 
weak-coupling (BCS) superfluidity to Bose-Einstein con- 
densation (BEC) of bound diatomic molecules. The key 
to this development is that, due to Feshbach molecular 
resonances, the strength and sign of the effective inter- 
atomic interaction can be tuned by varying the exter- 
nal magnetic field. For positive scattering lengths the 
resonantly enhanced interaction has been used to cre- 
ate a long-lived gas of diatomic molecules and Bose- 
Einstein condensates of molecules. For a negative scat- 
tering length, experimental evidence has been obtained 
for an atom-pair correlated state analogous to the BCS 
superconducting state at weak coupling |jj.y|. Additional 
studies of collective modes @, |j| and the spectroscopic 
pairing gap y| in the crossover region also display be- 
havior consistent with pairing, the most convincing ex- 
periment to date being the observation of a vortex lattice 
in the entire BCS-BEC crossover region jfj. 

The simplest model describing the crossover is due 
to Eagles and Leggett 0, la, E| in which the effective 
fermion-fermion interaction is parametrized by the scat- 
tering length, a. This approximation for the interaction is 
valid as long as the effective fermion-fermion interaction 
does not vary significantly on the scale of the Fermi en- 
ergy, Ey , the energy relevant for the many-body problem 
jlClL flll| . In practice, large scattering lengths are realized 
by using Feshbach resonances, and resonances for which 
this condition holds are referred to as "broad". For such 
resonances, the atom- molecule coupling and the width 
parameter for the resonance are large. In the opposite 
case, the resonance is referred to as "narrow". Most ex- 
periments on the crossover to date have been performed 
with broad resonances. For the resonance in 6 Li, the ex- 
)erimental data on the evolution of condensate profile 
12l agrees well with the universal single-channel model 
13l Il4| . In 0, |j] measurements of the collective axial 



and radial modes of a trapped 6 Li gas were reported to 
be quantitatively in agreement with results obtained from 
the zero temperature BCS-BEC crossover (lflLllfl|. 

In the future, one may anticipate that experiments will 
be made on the crossover for narrow resonances. It is 
therefore of interest to investigate the influence of en- 
ergy dependence of the interaction on the crossover. The 
leading contributions to the effective interaction beyond 
the scattering length are expressed in terms of the effec- 
tive range, r e , which is defined in terms of the s-wave 
scattering phase shift 5 by the equation 

0(fc 4 ), 
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where k is the wave number for the relative motion. The 
influence of the effective range is also interesting from 
the point of view of the BCS-BEC crossover produced 
by varying the density of the gas [13. In this case, the 
effective interaction between atoms at the Fermi surface 
depends on the Fermi momentum, and hence depends 
on density. Several approaches to study effective range 
dependence are based on microscopic interaction poten- 
tials, e.g., coupled square wells []3], the modified Poschl- 
Teller potential [19], and the Gaussian potential [l7ll2(l|. 
For monotonic attractive potentials the effective range is 
positive and can be shown to be related to an average 
range of the microscopic potential. The effective range 
expansion applied to the multichannel scatterin g p rob- 
lem yields effective ranges of either sign. In Refs. |2lLl22| 
the properties of a square well plus square barrier model 
potential were matched with the scattering properties of 
a renormalized two-channel resonance model [2]| and the 
zero temperature atom-molecule BCS-BEC crossover was 
re-expressed in terms of an effective single-channel po- 
tential model with a large and negative effective range 

The purpose of this paper is to investigate the depen- 
dence of the crossover on the energy dependence of the 
effective two-body interaction. We shall start from a mi- 
croscopic model, a square-well, square-barrier potential, 



which has previously been employed in Rcf. [21]. This 
can be used to construct effective interactions with dif- 
ferent scattering lengths and an effective range that can 
be either positive or negative. We then solve the Eagles- 
Leggett equations, determining the gap and the chemical 
potential self-consistently. Non-trivial corrections due to 
the negative effective range have not been fully investi- 
gated so far and such calculations provide a useful guide 
to explore the BCS-BEC crossover for narrow resonances 
in the future. 

The paper is organized as follows: Section|n]contains a 
general introduction to the low energy effective range ex- 
pansion of the effective interatomic interaction in a dilute 
Fermi gas |2J,[2ll. In Section ITTTl the two-body proper- 
ties of the square-well, square-barrier potential are de- 
rived. Subsequently, in Section HVl we numerically solve 
the zero temperature BCS-BEC crossover problem using 
the Eagles-Leggett variational state 0, la] and re-express 
the crossover properties in terms of the low energy two- 
body properties. 



II. EFFECTIVE RANGE EXPANSION 

Historically, the effective range expansion originated in 
attempts to understand nuclear forces from low energy 
nucleon scattering experiments [24L l2q. Classical prob- 
lems such as nucleon-nucleon scattering, the deuteron 
spectrum, and a decay were studied by such methods. 
Recently, the effective range was introduced to the field 
of ultracold atom gases where its size and sign are related 
to the width of the Feshbach resonances. 

The dominant scattering process in an ultracold two- 
component Fermi gas is the binary elastic s-wave scat- 
tering between the atoms. When the typical wavelength 
A corresponding to the collision energy E ~ h 2 /mX 2 (m 
is the mass of the fermion atom) is much larger than 
the range of the inter-atomic potential U(r), this pro- 
cess, in the center-of-mass frame, can be described by 
the Schrodinger equation with a radial wave function rip 
of the s-state: 



-V 2 0(r) + t/(r)0(r) = Eip(r). 



(2) 



Let us consider an interatomic potential with a finite 
range r\, i.e., U(r) = at r > r\. Using the solutions 
■01 and -02 of Eq. (2J with the energies Ei = h 2 k 2 jvn 
(i = 1,2), we get 
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K) / 0i0 2 dr. (3) 



Following Bethe [2jj, we introduce a reference wave func- 
tion ip Te ii which coincides with ipi outside the potential: 



ipi 



sin(fc i r + (5 l ), 



(4) 



where Si is the phase shift for energy Ei . For the reference 
wave functions, a similar relation to Eq. <J2j holds 
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Subtracting Eq. © from Eq. (0 leads to 
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Then we take fei = and use the relation 

1 
a 



lim k\ cot 8\ 

fci^O 



for the s-wave scattering length a to write 



kcotS = h -k 2 p(k), 

a 2 



with 



p(k) = 2tdr(^^ 

J \ sin o sin o 
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(8) 
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here we use the subscripts zero for the state with k = 
and we omit the subscripts 2. 

Comparing Eqs. (JJ and ((HJ, we get 



r e = p(k = 0) = 2 / dr 1 



n 



■ 2 r 

sm do 



(10) 



On resonance (sin<5o = 1), the above equation can be 
written as 



r e = r e0 = 2 / dr(l -0 O ), 
Jo 



(11) 



here we introduce r e o denoting the effective range at res- 
onance. 

In the case of the square-well potential with depth 
Uq = h 2 ko/m, ipo(r) = sin(fcor) for r < r\ with 
^o r i = (n + h)ir for (n — 0, 1, 2, • ■ ■ ) at resonance. Thus 
we get r e o = n for the square-well potential. However, 
the right hand side of Eq. il 111 is different from n for 
general potentials with a finite range. As will be seen 
later, both positive and negative r e o can be obtained us- 
ing a square-well, square-barrier potential. 

On resonance (So = n/2), the amplitude of the wave 
function for r < r\ becomes significant and it can be 
comparable to or greater than that outside the potential. 
From Eq. Ijllll we see that, while for positive r e o the am- 
plitude inside the potential is suppressed relative to that 
outside, for negative r e o the former is enhanced compared 
to the latter. Although the effect of r e is small away from 
resonance due to the first term on the right hand side of 
Eq. QJ being much larger than the second one, the sign 
of r e has important implications for the behavior of the 
wave function at resonance. 




Figure 1: Schematic picture of the square- well, square-barrier 
potential. We take Uo,Ui > in this paper. 



III. SQUARE- WELL, SQUARE-BARRIER 
MODEL POTENTIAL 



We shall model the microscopic interaction by a single- 
channel potential with a square well and a square bar- 
rier. The advantage of this potential, is that it contains 
the physics necessary to describe resonances: by increas- 
ing the thickness of the barrier, or its height, tunneling 
through the barrier is decreased, and therefore the po- 
tential can model the effects of changing the strength of 
coupling between continuum states and resonant states 
in the attractive well. 



A. Scattering length and effective range 

Let us consider two atoms interacting through a step- 
wise potential with an attractive core of radius r$ > 
and strength — Uq < and which in addition has a re- 
pulsive barrier of width w — r\ — ro for rx > r® and a 
strength U\ > (see Fig. QJ: 



-U for < r < r , 
U(r) = < U\ for r < r < r 1 , 
for n < r < oo . 



(12) 



This problem is solved by the standard method of match- 
ing the wave function solutions and their radial deriva- 
tives for the above three regions at their respective 
boundaries, rp, and r\ J2fj. 



The general solutions for the three regions are 

{Asin(Kr) for < r < ro, 

B (sinh(Kr) + Ccosh(«r)) for ro < r < n, 
sin(fcr + S) for n < r < oo. 

(13) 
The well depth parameter is conventionally defined as 
Kq = mUo/h 2 , the barrier height parameter as k\ = 
mlli/h 2 , and the relative kinetic energy of the two collid- 
ing particles is k 2 = mE/ti 2 . Furthermore, we introduce 
the auxiliary parameters K 2 = m(Uo + E)/h 2 = k% + k 2 , 
and k 2 = (Ui—E)/h 2 = k\ — k 2 . The energy dependent s- 
wave scattering phase shift is denoted by S(k). Equating 
the logarithmic derivatives at the boundaries, the follow- 
ing solutions are found for C and 5. First, C is found 
from the boundary condition at tq to be 



C 



Ktan(-KYo) — lTtaiih.(/wo) 

K — k tan(ii>o) tanh(«ro) 



(14) 



and the scattering phase shift is found from the boundary 
condition at r± 



' k tanh(wi) + C 
o(k) = — kri + arctan — 



= — kr\ 



arctan(7£), 



Ctanh(«ri) 



where 



TZ{k) = 



k Ktan(lvro) + K tanh[K(ri — ro)] 



K K + Ktan(AVo) tanh[«(ri — ro)] 
The scattering length is given by 
tan<5(/c) 



(15) 



(16) 



a = — lim 

fc^O 

1 



k 



=ri - - r — {k\ tan(fc r ) 



fc tanh[fci(ri -ro)]}. 



where 



( = k + kx tan(fc r )tanh[/si(ri — ro)]. 



(17) 



(18) 



In the present work we are also interested in the effective 
range, which can be derived from the low energy expan- 
sion QJ of the phase shift. Expanding kcotS for Eq. 
ltl5)l to the first order in energy, the first term yields the 
inverse scattering length as calculated above, and from 
the second term the following expression for the effective 
range is obtained: 



rx 



(19) 



with 



_ _ fcp +fci ^ 
k Q k 2 (a() 



^. S ech 2 [kx(rx-r )] 



= kl t kl A il - J-tanh[fcx(n - r )] 



kok 2 (aQ a \ kxrx 
1 r? 



k 2 a 3a 2 ' 



(21) 



At resonance where a diverges (the combination of a£ 
stays finite), the vanishing part r v is zero, and the ex- 
pression for r e reduces to 



>Vo =r\ 



k Q 



k\ 



tanh[fci(?'i -To)] 



k 2 fci 1 - tanh 2 [hi (n - r )] 
x {1 + fcir tanh[fei(n - r )]} , 



or equivalently 



r e o = n. + 



fc 



A 2 



tan(fc r ) - fc r 



fco 



fc 2 tan 2 (fcoro) — fc 2 , 



(22) 



(23) 



In addition, the condition for resonance (C = 0) becomes 
tan(fcoro) < —ko/ki. We can see, from Eq. (1221 . that the 
effective range r e o at resonance is smaller than T\ since 
< tanh[fci(ri — ro)] < 1. In the weak barrier limit, 
£ = fci(ri — ro) <C 1, the effective range at resonance can 
be written as 



r e o^r --i£ + 0(£ 2 ). 



(24) 



Thus adding a weak barrier to a square-well potential 
with the range r yields smaller effective range at res- 
onance than that for the original square-well potential 
r e0 = tq . The decrease of r e o means an enhancement of 
the amplitude in the potential region, which can be un- 
derstood as an effect of blocking probability leakage due 
to the additional barrier. 



B. Effective range and scattering properties 

To begin with, we illustrate the physical meaning of 
the effective range explicitly. In Fig. we plot the wave 
functions for \a\ —* oo in the cases of (a) positive and 
(b) negative effective range for the same values of k, ro, 
ki, and different r\. It is clearly seen that the resonance 
wave function for a negative effective range has a large 
amplitude in the potential range as mentioned before. 
Even in the case of Fig. El (a), a small bump can be 
observed close to the edge of the potential, which means 
a small decrease of the effective range compared to the 
square- well potential shown by Eq. <24ll . At r e = 0, the 
mean value of \^\ 2 for r < n is unity, rj -1 fj" 1 \i/j\ 2 dr = 1, 
thus / Q ri (l — ip 2 )dr = as given by Eq. (TTTTt . 

We now discuss the scattering properties of the square- 
well, square-barrier potential in detail. In Figs. and 
BJ we plot the parameter regions in which the effective 
range is negative and the absolute value of the scattering 
length is large close to resonance. Figure shows the 
dependence of these regions on the barrier width r\ — ro 
at fixed ro and fci. As n — ro increases, the width of 
the dark gray regions becomes smaller, which means the 
resonance becomes narrower. The effective range at res- 
onance decreases from r , which is for the square- well 
potential, to — oo with increasing r\ — r . This can also 
be seen from Eq. 12211 where the absolute value of the 
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Figure 2: Resonance wave functions for (a) positive and (b) 
negative effective range. In the both cases, we take kro = 0.01 
and fciro = 2.0, and ip is normalized so that V , ( r ) = sin(fcr + <5) 
for r > ri. 



second term increases monotonically from zero to infinity 
as r\ — ro increases with ro, fco, and k\ fixed. Figure 21 
shows the dependence of the above regions on the bar- 
rier height k\ at fixed ro and r\. Similarly to the above 
case, as k\ increases, the resonance becomes narrower and 
r e at resonance decreases from ro to — oo monotonically. 
This continuous decrease can be proved analytically by 
calculating the derivative of Eq. IJ23J) with respect to 
k\ and by taking the necessary condition for resonance, 
tan(fcoro) < —ko/ki, into account. From both Figs. 
and 01 we see that the resonance for the deeper potential 
well (larger ko) is narrower. 

Let us finally mention the behavior of the scattering 
length and the effective range away from resonance. In 
Figs. 0and|Hl we plot these quantities as functions of the 
width n — r and the height parameter ki of the barrier, 
respectively. The scattering length, which almost corre- 
sponds to the background one, increases linearly with r± 
for large barrier widths (see Fig. EJ: 



1 



a~ ri (25) 

However, this quantity does not change so much when 
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Figure 3: Parameter regions of fco and n — ro at fixed fciTo = 2 
for negative effective range (light gray region) and for reso- 
nance (dark gray region). Here the resonance region is defined 
by |ro/a| < 0.3; ro/o = —0.3 on the left boundary of the dark 
gray area and ro/a — 0.3 on that of the right one. The line 
plotted in the dark gray region corresponds to l/|o| = 0. The 
boundary of the light gray region corresponds to r e = 0, and 
the curve in this region shows r e = — oo. 



fei increases. It starts from a value comparable to r% at 
k\ = and converges to n as fci — > oo (see Fig. UJ|. 
Whether it approaches r± from above or below depends 
on the value of fcoHj. The behavior of the effective range 
is similar to that of the scattering length. It increases 
linearly with n at large n — ro (Fig. EJ : 



Figure 4: The same as Fig. EJfor parameter regions of fco and 
fci at fixed r 1 — 2r for negative effective range (light gray 
region) and for resonance (dark gray region). 
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(26) 
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and, with increasing fci, it remains of order n and con- 
verges to 2n/3 as fci — + oo (Fig. EJ. 



Figure 5: (Color online) The scattering length a (solid line) 
and the effective range r e (dash-dotted line) as functions of 
the barrier width n — ro. We take koro — 2 and fciro = 2. 



Basic formalism 



IV. NON-UNIVERSAL ZERO TEMPERATURE 
BCS-BEC CROSSOVER 



In the present section we examine the BCS-BEC 
crossover problem for a gas of equal numbers N/2 in two 
different internal states interacting through the square- 
well, square-barrier potential within an effective single- 
channel model. 



The Hamiltonian H for the two-component Fermi gas 
H = Y^ £kai a a k „ 

rj 2J ^kk'aj^alkja-k/jakT, (27) 



k.rr 



kk' 



,t 



where 0,^,(1^ are the creation and annihilation opera- 
tors, respectively, of atoms with wave vector k, and pseu- 
dospin index, a = | or J. . The single particle kinetic 




Figure 6: (Color online) The scattering length a (solid line) 
and the effective range r e (dashed line) as functions of the 
barrier height fei . The dotted lines show the asymptotic values 
of these quantities. We take feoro = 2 and ri = 2ro. 



energy is £k = ft 2 k 2 /(2m), and the inter-atomic interac- 
tion in momentum space is £/kk'- The zero temperature 
BCS-BEC crossover problem is here numerically studied 
along the lines of Refs. @, ED, El HI E3 based on the 
BCS variational state 



*bcs> = n k (u k + fka kT al k J |*o 



(28) 



where |*o) is the vacuum state, and u^ and Uk 
are the variational parameters normalized according 
to the condition |uk| 2 + |«k| 2 = 1. The energy 
of the BCS state E(iik,Vk) = (*bcs|^|*bcs) is to 
be minimized with respect to the variational param- 
eters Uk and «k, subject to the normalization condi- 
tion and the requirement that the mean particle num- 
ber N(u k , wt) = Ek,a(*BCs|aL a ka|*BCs) = 2 £ k \v k \ 2 
should be fixed. By introducing the Lagrange multipliers 
i?k and /j, to ensure normalization and number conserva- 
tion, the stationary condition 8[E — fiN — J2k ^k(|uk| 2 + 
|u k | 2 )]/<5(wk,Wk) = yields £ k it k + A k «k = ^k^k and 
MkA£ - £ k "k = E k Vk, where £k = h 2 k 2 /(2m) - \x. 
Here, the momentum dependent gap function is defined 
as Ak = J2k' ^kk'Wk'W k /- By solving the above equa- 
tions, one determines the usual quasi-particle energy 
Ek = (£k + A 2 .) 1 / 2 and the coherence factors |u k | 2 = 
(1 + £k/£k)/2 and \v k \ 2 = (1 - gk/^k)/2. Inserting the 
expression for u k v k — Ak/(2£k) into the definition of 
the gap yields the gap equation: 



1 v^ 



Al 



2E U , 



(29) 



This is to be solved self-consistently together with the 
equation for the atomic number density 



y?0 



in 

E k 



(30) 



At fixed density n = fcp/(37r 2 ), where k^ is the Fermi mo- 
mentum for the non-interacting gas, Eq. IJ30J1 determines 
the chemical potential, /U. For the contact interaction, 
the regularized homogeneous BCS-BEC crossover prob- 
lem [3, La| is analytically solvable [28L l29L I3(l | and results 
in the zero range (Leggett) reference curves in Figs. 00 

[EUD21 andUJ 

The partial wave decomposition of the interaction po- 
tential is 

Uw = J d 3 r e^'^r U(r) 

= ^(2^+l)47^C/ £ (fc,A:')P£(fc■fc , ), (31) 



e=o 



with 



U t (k,k')= drr 2 U {kr)U{r)j t {k'r). (32) 



Here jt(x) = y/Tv/{2x)j£+i/2(x) is the ^'th spherical 
Bessel function defined from the i/'th ordinary Bessel 
function of the first kind, J v (x), and Pi is the £'th Legen- 
dre function with k-k' = cos 9, 9 being the angle between 
the unit vectors k and k . We also used the addition for- 
mula {21 + l)Pt(k ■ k') = ^Y? m =^Yl m (k)Y tm {k') for 

the spherical harmonic functions Yt m {k). The momen- 
tum dependent gap can be decomposed in the same way 



A k = ^2 5Z Atrnjl(k)Yt m (k), 



(33) 



1=0 m^ 



where the A^ m 's are the angular weight coefficients. With 
the assumption that the main low energy contribution to 
the gap is of s-wave (£ = 0) character, the gap equation 
reduces to 



A(fc) = - 

IT 



dh , k , 2 U Q (k,k')A(k') 



where the s-wave part of the potential is 

/>OC 

U (k, k') = / dr r 2 jQ (kr)U(r)j (k'r) 
Jo 



(34) 



Uo + Ui 
2kk' 



sin(|fc + fc'|r ) sm(\k-k'\r ) 



2kk' 



\k + k'\ \k-k'\ 

sin(\k + k'\ri) sinflfc - k'\n) 



\k + k> 



\k-k' 



(35) 



B. Solution of the crossover equations 



The crossover equations consist of Eq. IJ30J1 and Eq. 
Il34t together with the expression (I35|l for the momen- 
tum dependent interaction. We notice that, for a finite 



range interaction potential vanishing beyond a scale 7*1, 
all momentum integrals are cutoff typically at the scale of 
1/Yi and therefore remain ultraviolet convergent. Here, 
r\ characterizes the range of the interatomic interactions. 
In the context of ultracold atoms, r\ can be taken to be 
the characteristic scale of the van der Waals interaction. 
We shall consider a system at fixed density and will study 
the properties of the system as the potential is changed by 
varying the well depth, the barrier height or the barrier 
width. Finally, in order to describe the crossover in terms 
of physically relevant two-body quantities, we paramet- 
rically plot quantities versus the scattering length a ob- 
tained from Eq. IJ17II and the effective range r e obtained 
from Eq. ltl9|l as functions of the model parameters. 

Below, we briefly describe our numerical procedure 
used to solve the crossover equations. By discretizing 
the momenta on a finite grid, the gap equation IJ34T) and 
the number density equation ll-'iOt were solved iteratively 
for the given density. The spacing of the momentum grid, 
Afc, was chosen to provide sufficient sampling on the scale 
of fcp ~ n 1 ' 3 . The maximum momentum fc max was cho- 
sen to be much larger than the scales of the internal range 
of the potential r^ 1 and 7*7/ . All results shown below are 
calculated with Afc = fcp/15 and fc ma x = 2500Afc (this 



corresponds to fc n 



50r for the model parameters 



employed). The system was driven from the BEC to the 
BCS limit and the initial guess for the fc-dependent gap 
function A(fc) was determined from the Fourier transform 
of the real space wave function determined in Section lTTTI 
The equations were solved by iteration, and the process 
was continued until the difference in /j, between succes- 
sive iterations was less than one part in 10 6 . The results 
were checked to be insensitive to further increase of grid 
size, fcmax, and decrease of grid spacing, Afc. 

We consider the case of driving the BCS-BEC crossover 
by varying the potential well parameter fco around its 
critical value fco c at which the scattering length di- 
verges. Both the scattering length and the effective 
range depend on fco and therefore r e /ro depends on 
(fcFCi) -1 through the parameter fco. Figure contains 
parametric plots of the effective range versus the scat- 
tering length for different values of the barrier width 
7*1 — r*o. As we are interested in the scattering prop- 
erties of the model and not in the particular values 
of the model parameters, it is more natural to show 
the r e curves in terms of the effective range r e o at 
(^Ffl) -1 = 0. The barrier widths used in the calculations, 
(7*1 -7*o)/r = 0.1,0.25,0.30,0.35, and 0.40, correspond 
to effective ranges r e o/ro sa 0.78, 0.12, —0.23, —0.66, and 
-1.20. The remaining parameter fci was chosen to be 
fciro = 2.0, and the density was chosen as n — 0.001r ( 



-3 

> 



i.e., fcpro — 0.31. For relative large and positive r e o, the 
effective range is only weakly dependent on the coupling 
constant in the crossover region as evidenced by the solid 
line in Fig. As r e o decreases, the effective range devel- 
ops a strong dependence on the inverse coupling constant. 
We note that similar results for r e presented in Fig. Q 
could have been obtained by changing the barrier height 
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Figure 7: (Color online) The effective range r e /r as a func- 
tion of the dimensionless constant (fcFa) -1 . The density is 
set to be fcpro ~ 0.31. 



instead. Such a degeneracy is to be expected for poten- 
tials with many model parameters since the mapping of 
the model parameters to the effective low energy param- 
eters is a many-to-two mapping. In the Eagles-Leggett 
model, only the scattering length appears and it is there- 
fore interpreted as the r e o = case with r e = for all 
couplings. 



C. Energy gap 

We first examine the effective range dependence of the 
gap parameter A(fc). In panels (a)-(c) of Fig. HJ A(k)/E-p 
obtained from the crossover equations is displayed for 
several values of r e o, here E-p = ft 2 fcp/2m is the Fermi 
energy for the non-interacting gas. The solid curves in 
panels (a)-(c) show the evolution of A(fc) in the near BCS 
limit. In this limit, the gap has a maximum at fc = for 
relatively large and positive r e0 as displayed in panel (a) . 
As r e o decreases, the maximum flattens and the sign of 
the curvature at fc = changes from negative to positive 
as can be seen by comparing the solid curves in panels (a) 
and (b). For negative r e o, the gap has a local minimum 
at fc = and the maximum of the gap moves to a higher 
fc [see the solid curve in panel (c)]. 

The energy gap in the ground state is the minimum 
of the quasiparticle dispersion Ek for all k yj and, for 
a Fermi gas with contact interactions, the gap Ak = A 
is momentum independent. For a positive /j,, the energy 
Ek = (£k + A 2 ) 1 / 2 is minimized at h 2 k 2 /2m = fx with 
the value E min = |A|. On the other hand, for a negative 
/i, the energy minimum -Emm = (/-t 2 + A 2 ) 1 ' 2 is attained 
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Figure 9: (Color online) The chemical potential h/Ef ver- 
sus (kFa)~ is plotted for changing the depth of the potential 
well j22|. Each curve is labeled by the value of the dimen- 
sionless effective range r e o/ro in the unitarity limit. As the 
effective range becomes negative, the finite range correction 
to the universal chemical potential changes its sign. In the in- 
set, we have plotted the chemical potential fio in the unitarity 
limit, (fefl) - = 0, as a function of the unitary effective range 
r e o. The universal result /xo/-Ef ~ 0.59 obtained by Leggett 
Q is marked with the horizontal dashed line. The density is 
set to be k F ro ~ 0.31. 



Figure 8: (Color online) The solution of the gap equation for 
the values (fed) -1 = —1,0, and 1 of the inverse coupling 
constant is presented for several values of the effective range 
r e o in the unitarity limit. For positive effective range (a), 
the zero momentum gap is a local maximum. Close to zero 
effective range, the maximum turns into a plateau with a 
vanishing momentum gradient. In (c), the zero momentum 
gap is turned into a local minimum for all values of the inverse 
coupling. The density is set to be k F ro — 0.31. 



Next we consider the functional behavior in the gen- 
eral case of a momentum dependent gap Ak- We assume 
that the Fermi gas is so dilute that the range of the po- 
tential is much shorter than the interatomic spacing, i.e., 
r o, r i "C fcp . The momentum dependent gap varies on 
the scale of 1/Vo which is large compared with kp, and 
therefore the behaviors described above for a constant 
gap are expected to be rather general |8lT |. 



at k = 0. Therefore, as the interaction is decreased from 
the BEC limit to the BCS one, the fc m j n value, where E^ 
has its minimum, shifts from to y / 2mfi/h 2 at \x = 
to fc m in = k-F in the BCS limit. Due to the interaction- 
induced shift in the position of k m i n , the functional be- 
havior of Ek is also directly reflected in the pair wave 
function ip(k) = A/(2_Ek) which appears directly in the 
gap equation, and also in a number of derived quantities 
such as the healing length and condensate fraction. In 
the BEC limit, where the gap equation reduces to the 
equation for two bound particles, i/)(k) represents the in- 
ternal wave function of the pair. The l/E k dependence 
implies a maximum of ^(fc) at k = for /i < (the BEC 
case) and at ^/2rnjl/h for \i > 0, the BCS case (see also 



D. Chemical potential 

Next we discuss the influence of the nonzero effective 
range on the chemical potential. Figure|§| contains a plot 
of the chemical potential as a function of (kpa)^ 1 for sev- 
eral values of r e o . At a sufficiently low value of r e o , the 
chemical potential approaches the universal result (short 
blue dashes denoted by "Leggett") as obtained from the 
zero-range model. By an additional variation of the bar- 
rier parameter, the effective range r e o/ro is varied from 
a positive value 0.78 (solid line) to —1.20 (dotted line). 
The plots of the chemical potential for different values 
of r e0 are all seen to cross at (fcpa) -1 ~ 0.55, where n 
is close to zero. The inset in Fig. 0J contains a plot of 
the chemical potential jjlq at fcp|a| ^ °o as a function 



of the effective range, which clearly shows the correction 
due to the effective range. The horizontal dashed line of 
Ho/Ep » 0.59 corresponds to jlto for a contact potential, 
and the crossing of this line at r e $ ~ also validates our 
approach to discuss non-universal corrections by use of a 
finite-range model potential. For positive r e0 , the chem- 
ical potential for kp\a\ — > oo is enhanced with respect 
to the zero range case as was also observed in Ref. J2fJ 
based on the Gaussian potential. On the other hand for a 
negative r e0 the chemical potential is reduced compared 
to the curve for a contact potential. These non-universal 
corrections become small with decreasing density. We 
have observed that, for a lower density of n = O.OOOlr^" , 
i.e., kpro ~ 0.14, deviations from the Leggett curve in 
the region —1 < (kpa)" 1 < 1 are about a factor of two 
smaller than those in Fig. 

By matching the scattering properties of single- and 
two-channel models on resonance, it is found that 
4nH 2 a/m = —g 2 jv and r e0 = — %Tth A /(mg 2 ), where g is 
the atom-molecule coupling and v is the detuning in the 
two-channel boson-fermion model [231. Because of r e0 is 
negative in this model, for (kpa)^ 1 ~ 0, the above non- 
universal reduction of the chemical potential is expected 
in less wide to narrow Feshbach resonances [2ll l22|. It 
is also noted that this reduction close to a narrower res- 
onance cannot be achieved with a monotonic attractive 
potential (which all have r e Q > 0) and therefore con- 
strains the model potential when modelling atom gases. 

In Fig. E3 we compare the chemical potential and 
the binding energy E\, of the two-body bound state 
of the same square-well, square-barrier potential. For 
a ~^> ro,T*i, the binding energy is well approximated by 
using only a and r e as Eh ~ (H 2 /ma 2 )(l + r e /a) (Eh is 
defined to be positive), and thus the negative effective 
range reduces the binding energy close to threshold. For 
a < 7*0) ?"i, since the extension of the wave function of the 
bound state is comparable to the range of the potential, 
Eh depends strongly on its specific shape. Therefore, Eh 
for a potential with a larger negative r e o, which has a 
stronger barrier [i.e., larger £ = fei(ri — ro)], is smaller at 
(fcFa) -1 3> 1- We can also see that the chemical potential 
consistently behaves as \x — ► —Eh/2 in this limit. 



E. Consequences of the nonzero effective range in 
non-universal resonances 



Most experiments on the BCS-BEC crossover done so 
far have been performed with broad resonances. Accord- 
ing to an estimate in J22|, kpr e ~ 0.01 and r e /ro < 1 (23 
in the experiment by Regal et al. y|, and kpr e ~ 10 -4 
and r e /ro ~ 0.1 in the experiment by Zwierlein et al. Q. 
Due to the small values of kpr e , the effective range has 
an imperceptible effect in these experiments. However, it 
is essential to study narrower resonances, in which effec- 
tive range corrections are significant, for understanding 
strongly interacting Fermi gases more deeply. Therefore, 
consequences of the nonzero effective range, especially 
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Figure 10: (Color online) Comparison between the chemi- 
cal potential (solid lines) and the binding energy Eh of the 
two-body bound state (dashed lines) for several values of the 
effective range r e o/ro = 0.78 (red), 0.12 (green), —0.23 (blue), 
and —1.2 (magenta). The density is set to be k?r<z ~ 0.31. 



for negative cases, in future experiments on narrow res- 
onances are worth discussing. In the remaining part of 
this paper, we study effective range corrections for uni- 
form systems on the following four measurable quantities. 
zero momentum gap, momentum distribution, and con- 
densate fraction. 



1. Zero momentum gap 

In Fig. ^]the zero momentum gap A(0) is plotted as a 
function of (kpa) -1 . All curves cross at almost the same 
point near (kpa)^ 1 = and thus the zero momentum gap 
shows little dependence on r e o close to resonance. In the 
BCS region [(fcpa) -1 < 0], however, the effective range 
corrections can be significant and A(0) is enhanced for 
large and negative r e o. A larger negative value of r e o, 
which is yielded by a stronger barrier [larger £ = k\ (r\ — 
ro)] in the present single-channel model, corresponds to a 
weaker atom-molecule coupling g in a two-channel model. 
Thus, in this case, BCS pairing is more favorable on the 
BCS side and the gap is enhanced. This result can be 
also understood in terms of the effective range expansion 
introducing the effective scattering length a e (k): 



-J- = 1 - l -r e k 2 . 
a e (k) a 2 



(36) 



For negative a and r e , the absolute value of the effective 
scattering length |a e | is larger than \a\; i.e., the effective 
range correction due to a negative r e yields a stronger 
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Figure 11: (Color online) The zero temperature gap A(0)/-Ef 
at k = is plotted as a function of the inverse coupling con- 
stant (fea) -1 for several values of the effective range r e o/ro. 
The density is set to be k^ro ~ 0.31. 



attractive interaction on the BCS side, which leads to a 
larger BCS gap. Using the effective scattering length, 



we may write A(0)/E F ps (8/e 2 )e 



-2{/£ F a c (fe F )}" 1 / 7r - 

(8/e 2 )e{ 2 ( fcFa ^ 1 - fcFr ">/ 7r in the BCS limit, which con- 
verges toward (8/e 2 )e- 2(fcFarl/lr in the far BCS limit. 

In the far BEC regime [(kpa)^ 1 ^> 1], where the 
molecules are deeply bound, A(0) also strongly depends 
on the specific properties of the potential, as has been 
seen for the chemical potential, and therefore it deviates 
from the zero range result. The suppression of A(0) for 
a large, negative r e o compared to that for a small neg- 
ative or positive r e o can be attributed to the stronger 
repulsive barrier in the former case, which acts to de- 
stroy the pair in this regime. We have also confirmed 
that, with decreasing density, A(0) for each value of r e o 
approaches the one for a contact potential. For a lower 
density of n = O.OOOlr,^ , the effective range corrections 
in the region of —1 < (fcpa) -1 < 1 are about a factor of 
two smaller than those in Fig. El as in the case of the 
chemical potential [38]. 



2. Momentum distribution 



The momentum distribution rik/2 = |ffe| 2 /V of fermion 
atoms (for spin-up or spin-down) is shown in Fig. 
1121 where n^ being 2\vk\ 2 /V. Here we plot rifc/2 for 
(A^a) -1 = —1,0, and 1 for the same parameter sets as 
in Fig. 03 In the near to far BEC limits of (fcpa) > 1, 
all curves for different r e0 almost coincide with the one 
for the contact potential. However, from the unitarity 
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Figure 12: (Color online) The momentum distribution n*/2 
of fermion atoms for (fei) -1 = —1,0, and 1 is presented for 
several values of of the effective range r e o in the unitarity 
limit. The density is set to be kpro ~ 0.31. 



regime (fcpa) -1 ~ to the near BCS limit (fcpa) -1 ~ — 1, 
we observe significant effective range corrections. The 
results for (fcpa) -1 = —1 and in this figure show that 
corrections due to the negative effective range suppress 
nk=o and broaden the momentum distribution. As dis- 
cussed before, a large, negative r e o, which corresponds 
to a weak atom-molecule coupling, favors the BCS state 
on the BCS side. Therefore the Fermi surface is smeared 
and nk is broadened in this case. 



3. Condensate fraction 

We consider the influence of the effective range on the 
condensate fraction in the BCS-BEC crossover. In a 
Fermi gas, the appearance of superfluid order is related to 
the appearance of off-diagonal long range order as found 
in the two-body density matrix. Following [34l \WL l3q. 
we calculate the number of the condensed fermion pairs 
as N = X)kl u k^k| 2 and examine the influence of the 
nonzero effective range. 

FigurefL3l contains a plot of the normalized condensate 
fraction Nq/{N/2) = 2Nq/N of a Fermi gas as a function 
of the inverse coupling for the same parameter sets as 
those used in Figs. El and El Within the mean-field ap- 
proximation, the normalized molecular condensate frac- 
tion converges to unity in the BEC limit independent of 
the detailed nature of the potential. On the other hand, 
when the binding energy is decreased, the dimers dissoci- 
ate more easily and the condensate fraction correspond- 
ingly decreases and eventually vanishes in the BCS limit. 
For a large and negative value of r e o , the dissociation is 
suppressed due to the weak atom-molecule coupling, and 
thus the condensate fraction is enhanced as shown in this 
figure. 
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Figure 13: (Color online) The normalized zero temperature 
condensate fraction No/ (N/2) as a function of the inverse 
coupling constant (fed) -1 for several values of the effective 
range r e o/ro. The density is set to be kpro ~ 0.31. 



the appearance of negative effective range, which is pre- 
dicted by two-channel models. 

First, we studied the two-particle scattering problem 
with the square-well, square-barrier potential, and have 
clearly demonstrated the physical meaning of the nega- 
tive effective range. The effective range at resonance was 
carefully analyzed and the analytical expression for the 
weak barrier case was also presented. This expression 
shows that adding a barrier to a square-well potential 
yields smaller values of the effective range. 

We then applied this potential to the BCS-BEC 
crossover phenomena within an effective single-channel 
potential model. By solving the crossover equations self- 
consistently, we studied the effects of both the positive 
and negative effective range within a unified framework. 
We have observed that corrections due to the negative 
effective range, in general, appear as a maximum of the 
gap parameter at finite momentum and a reduction of 
the chemical potential close to the resonance and in the 
BCS region. In view of the possibility of performing ex- 
periments on systems with narrow resonances, we have 
discussed consequences of the nonzero effective range on 
various measurable quantities [23|. The results presented 
in this work will be helpful when one explores the BCS- 
BEC crossover for narrow resonances. 



SUMMARY AND CONCLUSION 



We have presented a detailed study of the two-body 
bound state and scattering properties of a simple solvable 
model potential with a barrier to model the narrow res- 
onance phenomena in the BCS-BEC crossover in atomic 
Fermi gases. We paid particular attention to the effect 
of the nonzero effective range and its sign as the inverse 
scattering length was tuned from large negative values to 
large positive ones. The advantage of our simple single- 
channel potential model is that it explicitly displays the 
relation between the potential resonance phenomena and 
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